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Abstract 



Yields of prompt and non-prompt ]/tp, as well as Y(1S) mesons, are measured by the 
CMS experiment via their u + }T decays in PbPb and pp collisions at y/s^ = 2.76 TeV 
for quarkonium rapidity \y\ < 2.4. Differential cross sections and nuclear modifica- 
tion factors are reported as functions of y and transverse momentum pj, as well as 
collision centrality. For prompt ]/ip with relatively high pj (6.5 < pj < 30GeV/c), a 
strong, centrality-dependent suppression is observed in PbPb collisions, compared to 
the yield in pp collisions scaled by the number of inelastic nucleon-nucleon collisions. 
In the same kinematic range, a suppression of non-prompt J/t/>, which is sensitive to 
the in-medium b-quark energy loss, is measured for the first time. Also the low-px 
Y(1S) mesons are suppressed in PbPb collisions. 
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1 Introduction 

At large energy densities and high temperatures, strongly interacting matter consists of a de- 
confined and chirally-symmetric system of quarks and gluons [1J. This state, often referred to 
as "quark-gluon plasma" (QGP) [2], constitutes the main object of the studies performed with 
relativistic heavy-ion collisions (IHSl. 

The formation of a QGP in high-energy nuclear collisions can be evidenced in a variety of ways. 
One of its most striking expected signatures is the suppression of quarkonium states [7j, both 
of the charmonium (J/tp, ip', Xc, etc.) and the bottomonium (Y(1S,2S,3S), Xb, etc.) families. 
This is thought to be a direct effect of deconfinement, when the binding potential between the 
constituents of a quarkonium state, a heavy quark and its antiquark, is screened by the colour 
charges of the surrounding light quarks and gluons. The suppression is predicted to occur 
above the critical temperature of the medium (T c ) and depends on the QQ binding energy. 
Since the Y(1S) is the most tightly bound state among all quarkonia, it is expected to be the one 
with the highest dissociation temperature. Examples of dissociation temperatures are given in 
Ref. II: T dissoc ~ 1 T c/ 1.2 T c , and 2 T c for the Y(3S), Y(2S), and Y(1S), respectively. Similarly, 
in the charmonium family the dissociation temperatures are < 1 T c and 1.2 T c for the xp' and 
]/xp, respectively. However, there are further possible changes to the quarkonium production 
in heavy-ion collisions. On the one hand, modifications to the parton distribution functions 
inside the nucleus (shadowing) and other cold-nuclear-matter effects can reduce the production 
of quarkonia without the presence of a QGP (9j[l0l. On the other hand, the large number of 
heavy quarks produced in heavy-ion collisions, in particular at the energies accessible by the 
Large Hadron Collider (LHC), could lead to an increased production of quarkonia via statistical 
recombination |[TTHT6ll . 

Charmonium studies in heavy-ion collisions have been carried out for 25 years, first at the 
Super Proton Synchrotron (SPS) by the NA38 HZI, NA50 HHCOO, and NA60 [H3 fixed-target 
experiments at 17.3-19.3 GeV centre-of-mass energy per nucleon pair (^/s NN ), and then at the 
Relativistic Heavy Ion Collider (RHIC) by the PHENIX experiment at v /s^ = 200 GeV EH. In 
all cases, ]/ip suppression was observed in the most central collisions. At the SPS, the suppres- 
sion of the ip' meson was also measured ||19 || . Experimentally, the suppression is quantified by 
the ratio of the yield measured in heavy-ion collisions and a reference. At RHIC, the reference 
was provided by the properly scaled yield measured in pp collisions. Such a ratio is called the 
nuclear modification factor, Raa- m the absence of modifications, one would expect Raa - 1 for 
hard processes, which scale with the number of inelastic nucleon-nucleon collisions. For bot- 
tomonia, the production cross section is too small at RHIC to make definitive statements [22 J. 
With the higher energy and luminosity available at the LHC, new studies for charmonia and 
bottomonia have become possible: (i) ATLAS has reported a suppression of inclusive ]/xp with 
high transverse momenta px in central PbPb collisions compared to peripheral collisions at 
■s/s^ = 2.76 TeV [23J; (ii) a suppression of the excited Y states with respect to the ground state 
has been observed in PbPb collisions at y/s^ = 2.76 TeV compared to pp collisions at the same 
centre-of-mass energy by the Compact Muon Solenoid (CMS) collaboration [24 J. 

At LHC energies, the inclusive ]/xp yield contains a significant non-prompt contribution from b- 
hadron decays |25l427| . Owing to the long lifetime of the b hadrons (0(500) f/m/ c), compared 
to the QGP lifetime (O(10) fm/c), this contribution should not suffer from colour screening, 
but instead may reflect the b-quark energy loss in the medium. Such energy loss would lead to 
a reduction of the b-hadron yield at high pj in PbPb collisions compared to the binary-collision- 
scaled pp yield. In heavy-ion collisions, only indirect measurements of this effect exist, through 
single electrons from semileptonic open heavy-flavour decays |28ti30| : to date, the contribu- 
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tions from charm and bottom have not been disentangled. The importance of an unambiguous 
measurement of open bottom flavour is driven by the lack of knowledge regarding key features 
of the dynamics of parton energy loss in the QGP, such as its colour-charge and parton-mass 
dependencies (31j|32l and the relative role of radiative and collisional energy loss 1133 1 . CMS is 
well equipped to perform direct measurements of b-hadron production in heavy-ion collisions 
by identifying non-prompt J/xp from b-hadron decays via the reconstruction of secondary 
vertices. 

The paper is organised as follows: the CMS detector is briefly described in Section [2] Section [3] 
presents the data collection, the PbPb event selection, the muon reconstruction and selection, 
and the Monte Carlo (MC) simulations. The methods employed for signal extraction are de- 
tailed in Section [4] Section [5] describes the acceptance correction factors and the estimation of 
the reconstruction efficiencies. The pp baseline measurements are summarized in Section [6] 
The results are presented in Section |7j followed by their discussion in Section [HJ 

2 The CMS Detector 

A detailed description of the CMS experiment can be found in Ref. |l34l| . The central feature of 
the CMS apparatus is a superconducting solenoid of 6 m internal diameter. Within the field vol- 
ume are the silicon tracker, the crystal electromagnetic calorimeter, and the brass/ scintillator 
hadron calorimeter. 

CMS uses a right-handed coordinate system, with the origin at the nominal interaction point, 
the x axis pointing to the centre of the LHC, the y axis pointing up (perpendicular to the LHC 
plane), and the z axis along the counterclockwise-beam direction. The polar angle 6 is mea- 
sured from the positive z axis and the azimuthal angle (p is measured in the x-y plane. The 
pseudorapidity is defined as t] = — ln[tan(0/2)]. 

Muons are detected in the interval \rj\ < 2.4 by gaseous detectors made of three technologies: 
drift tubes, cathode strip chambers, and resistive plate chambers, embedded in the steel return 
yoke. The silicon tracker is composed of pixel detectors (three barrel layers and two forward 
disks on either side of the detector, made of 66 million 100 x 150 ]im 2 pixels) followed by mi- 
crostrip detectors (ten barrel layers plus three inner disks and nine forward disks on either side 
of the detector, with strips of pitch between 80 and 180 /mi). The transverse momentum of 
muons matched to reconstructed tracks is measured with a resolution better than ~1.5% for pj 
smaller than 100 GeV/c II35I . The good resolution is the result of the 3.8 T magnetic field and the 
high granularity of the silicon tracker. 

In addition, CMS has extensive forward calorimetry including two steel/ quartz-fibre Cheren- 
kov forward hadron (HF) calorimeters, which cover the pseudorapidity range 2.9 < \rj\ < 
5.2. These detectors are used in the present analysis for the event selection and PbPb collision 
centrality determination, as described in the next section. Two beam scintillator counters (BSC) 
are installed on the inner side of the HF calorimeters for triggering and beam-halo rejection. 

3 Data Selection 
3.1 Event Selection 

Inelastic hadronic PbPb collisions are selected using information from the BSC and HF calorime- 
ters, in coincidence with a bunch crossing identified by the beam pick-up (one on each side of 
the interaction point) [34J . Events are further filtered offline by requiring a reconstructed pri- 
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mary vertex based on at least two tracks, and at least 3 towers on each HF with an energy 
deposit of more than 3 GeV per tower. These criteria reduce contributions from single-beam 
interactions with the environment (e.g. beam-gas collisions and collisions of the beam halo 
with the beam pipe), ultra-peripheral electromagnetic interactions, and cosmic-ray muons. A 
small fraction of the most peripheral PbPb collisions are not selected by these minimum-bias 
requirements, which accept (97 ± 3)% of the inelastic hadronic cross section II36I . A sample 
corresponding to 55.7 M minimum-bias events passes all these filters. Assuming an inelastic 
PbPb cross section of CpbPb = 7.65 b |36], this sample corresponds to an integrated luminosity 
of £i n t = 7.28 ub^ 1 . This value is only mentioned for illustration purposes; the final results are 
normalized to the number of minimum-bias events. 

The measurements reported here are based on dimuon events triggered by the Level-1 (LI) 
trigger, a hardware-based trigger that uses information from the muon detectors. The CMS 
detector is also equipped with a software-based high-level trigger (HLT). However, no further 
requirements at the HLT level have been applied to the LI muon objects used for this analysis. 

The event centrality distribution of minimum-bias events is compared to events selected by 
the double-muon trigger in Fig. [T] The centrality variable is defined as the fraction of the total 
cross section, starting at 0% for the most central collisions. This fraction is determined from 
the distribution of total energy measured in both HF calorimeters [37] . Using a Glauber-model 
calculation as described in Ref. f36l| , one can estimate variables related to the centrality, such 
as the number of nucleons participating in the collisions (N part ) and the nuclear overlap func- 
tion (Taa), which is equal to the number of elementary nucleon-nucleon (NN) binary collisions 
divided by the elementary NN cross section and can be interpreted as the NN equivalent in- 
tegrated luminosity per heavy ion collision, at a given centrality 1138 1 . The values of these 
variables are presented in Table [I] for the centrality bins used in this analysis. The double- 
muon-triggered events are more frequent in central collisions since the main physics processes 
that generate high-px muon pairs scale with the number of inelastic nucleon-nucleon collisions. 
In the following, N part will be the variable used to show the centrality dependence of the mea- 
surements. 




100 80 60 40 20 



Centrality (%) 

Figure 1: Centrality distribution of the minimum-bias sample (solid black line) overlaid with 
the double-muon triggered sample (hashed red) in bins of 2.5%. 
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Table 1: Average and root-mean-square (RMS) values of the number of participating nucleons 
(Npart) and of the nuclear overlap function (Taa) for the centrality bins used in this analysis [36 J. 

Npart Taa (mb _1 ) 



Centrality (%) Mean RMS Mean RMS 



0-10 


355.4 


8.3 


23.19 


1.12 


10-20 


261.4 


9.6 


14.48 


1.04 


20-30 


187.2 


7.3 


8.78 


0.72 


30^0 


130.0 


5.6 


5.09 


0.49 


40-50 


86.3 


4.3 


2.75 


0.31 


50-100 


22.1 


0.9 


0.47 


0.04 


0-20 


308.4 


6.3 


18.83 


0.77 


20-100 


64.2 


1.4 


2.37 


0.12 


0-100 


113.1 


1.7 


5.66 


0.18 



Simulated MC events are used to tune the muon selection criteria, to compute the acceptance 
and efficiency corrections, and to obtain templates of the decay length distribution of ]/tp from 



b-hadron decays. For the acceptance corrections described in Section 5.1 three separate MC 
samples without kinematic requirements are used: prompt ]/xp, ]/ip from b-hadron decays, and 
Y(1S). Prompt J/tp and Y(1S) are produced using PYTHIA 6.424 [39| at y/s = 2.76 TeV, which 
generates events based on the leading-order colour-singlet and colour-octet mechanisms, with 
non-relativistic quantum chromodynamics (QCD) matrix elements tuned BO! by comparison 
with CDF data [il]. The colour-octet states undergo a shower evolution. For the non-prompt 
]/ip studies, the b-hadron events are produced with PYTHIA in generic QCD 2— >2 processes. In 
all three samples, the J/i/> or Y(1S) decay is simulated using the evtgen [42] package. Prompt 
]/ip and Y(1S) are simulated assuming unpolarized production, while the non-prompt ]/xp po- 
larization is determined by the sum of the exclusive states generated by evtgen. Final-state 
bremsstrahlung is implemented using PHOTOS |43| . 



For some MC simulation studies, in particular the efficiency corrections described in Section 5.2 
each PYTHIA signal event is simulated with GEANT4 |44| and then embedded in a realistic 
heavy-ion background event. The background events are produced with the hydjet event 
generator [45 1 and then simulated with GEANT4 as well. The embedding is done at the level of 
detector hits and requires that the signal and background production vertices match. The em- 
bedded event is then processed through the trigger emulation and the full event reconstruction 
chain. Collision data are used to validate the efficiencies evaluated using MC simulations, as 
discussed in Section IBT21 



3.2 Muon Selection 

The muon offline reconstruction algorithm starts by reconstructing tracks in the muon detec- 
tors, called standalone muons. These tracks are then matched to tracks reconstructed in the 
silicon tracker by means of an algorithm optimized for the heavy-ion environment |46l |47| . 
The final muon objects, called global muons, result from a global fit of the standalone muon and 
tracker tracks. These are used to obtain the results presented in this paper. 

In Fig. [2} the single-muon reconstruction efficiency from MC simulations is presented as a func- 
tion of the muon pj and if 1 . The reconstruction efficiency is defined as the number of all re- 
constructed global muons divided by the number of generated muons in a given (jf , pj) bin. 
It takes into account detector resolution effects, i.e. reconstructed pj and tj values are used in 
the numerator and generated pj and t] values in the denominator. To obtain a clear separation 
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between acceptance and efficiency corrections, a detectable single-muon acceptance is defined 
in the (yf , pj) space. For the }/ip analysis this separation is defined by the contour that roughly 
matches a global muon reconstruction efficiency of 10%, indicated by the white lines superim- 
posed in Fig.|2j which are described by the conditions 



p\ > 3.4GeV/c 

p\ > (5.8-2.4 x |7f |)GeV/c 

p\ > (3.4 - 0.78 x |^|) GeV/c 



for \rf\ < 1.0, 
for 1.0 <\f] H \< 1-5, 
for 1.5 < |^| < 2.4. 



(1) 



Muons failing these conditions are accounted for in the acceptance corrections discussed in 



Section 5.1 Muons that pass this acceptance requirement can still fail to pass the trigger, track 
reconstruction, or muon selection requirements. These losses are accounted for by the efficiency 
corrections discussed in Section l52l 




Figure 2: Reconstruction efficiency of global muons in the (tj^ 1 , p^.) space, illustrating the lower 
limits (white lines) of what is considered a detectable single muon for the ]/xp analysis. 

For the Y(1S) analysis, where the signal-to-background ratio is less favourable than in the J/i/? 
mass range, a higher p l j is required than for the ]/xp analysis, 



p\ > 4GeV/c, 



(2) 



independent of t] f 



Various additional global muon selection criteria are studied in MC simulations. The MC dis- 
tributions of the ]/xp decay muons are in agreement with those from data to better than 2%, 
which is within the systematic uncertainty of the data/MC efficiency ratio (Section |5.2| . The 
transverse (longitudinal) distance of closest approach to the measured vertex is required to be 
less than 3 (15) cm. Tracks are only kept if they have 11 or more hits in the silicon tracker, and 
the x 2 per degree of freedom of the global (inner) track fit is less than 20 (4). The x 2 probability 
of the two tracks originating from a common vertex is required to be larger than 1%. From MC 
simulations we find that these criteria result in a 6.6%, 5.1%, and 3.9% loss of prompt ]/ip, non- 
prompt J/]/?, and Y(1S) events, respectively, not including the track reconstruction and trigger 
efficiencies. 
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4 Signal Extraction 

4.1 J/ip Analysis 
4.1 .1 Inclusive ]ltp 

The pair invariant-mass spectrum is shown in Fig. pun the region 2 < < 

4GeV/c 2 for muon pairs with < pj < 30GeV/c and rapidity \y\ < 2.4, after applying the 
single-muon quality requirements. No minimum pair-px requirement is applied explicitly. 
However, the CMS acceptance for pairs in this mass range requires a minimum pj that 

is strongly y-dependent and is ~ 6.5GeV/c at y = 0. The black curve in Fig. [3] represents an 
unbinned maximum likelihood fit to the p + }i~ pair spectrum, with the signal described by the 
sum of a Gaussian and a Crystal Ball function, with common mean mn and width a, and the 
background described by an exponential. The Crystal Ball function /cb(^) combines a Gaus- 
sian core and a power-law tail with an exponent n to account for energy loss due to final-state 
photon radiation, 



fcB(m) 
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exp 
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for 

f or m=m < 



> —a; 



-a. 



(3) 



The parameter a defines the transition between the Gaussian and the power-law functions. 
The fitted mean value, m = (3.090 ± 0.002) GeV/c 2 , is slightly below the PDG value of m ]/tp = 
3.097 GeV/c 2 [48J because of momentum scale effects in the heavy-ion data; the width is a = 
(39 ± 2) MeV/c 2 , consistent with MC expectations. The number of inclusive ]/ip mesons ob- 
tained by the fit is 734 ± 54. 
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Figure 3: Invariant-mass spectrum of Y- + }i~ pairs (black circles) with \y\ < 2.4 and 
< pj < 30GeV/c integrated over centrality. The fit to the data with the functions discussed 
in the text is shown as the black line. The dashed blue line shows the fitted background contri- 
bution. 



The analysis is performed in bins of the ]/ip meson pj and y, as well as in bins of event centrality. 
Integrating over all centrality (0-100%) and pi (6.5 < pj < 30GeV/c) the rapidity bins are 

\y\ < 1.2, 1.2 < \y\ < 1.6, and 1.6 < \y\ < 2.4. 
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For the two forward bins, the CMS acceptance extends to lower pj, so results are also presented 
for the bins 

1.2 < |y| < 1.6 and 5.5 < p T < 30GeV/c, as well as 1.6 < |y| < 2.4 and 3 < p T < 30GeV/c. 

These values allow a better comparison with the low-px measurements of the ALICE experi- 
ment, which has acceptance for ]/ip with pj > OGeV/c for the rapidity intervals |y| < 0.9 and 
2.4 < y < 4.0, in the electron and muon decay channels, respectively |49|. 

Integrating over all centrality (0-100%) and rapidity (|y| < 2.4) the pj bins are 

6.5 < p T < 10GeV/c and 10 < p T < 30GeV/c. 

Integrating over the pj range 6.5 < pj < 30GeV/c and rapidity \y\ < 2.4, the centrality bins 
are: 0-10%, 10-20%, 20-30%, 30-40%, 40-50%, and 50-100%. 

The unbinned maximum likelihood fit with the sum of Crystal Ball and Gaussian functions 
is performed in each of these bins. Because of the small sample size, the parameters of the 
signal shape are determined for each rapidity and pj interval, integrated over centrality, as 
the dominant effect on the mass shape is the pj- and rapidity-dependent mass resolution. The 
values are then fixed for the finer centrality bins. The background shape is allowed to vary in 
each bin. 

4.1 .2 Prompt and Non-prompt ]lxp 

The identification of ]/xp mesons coming from b-hadron decays relies on the measurement of 
a secondary p + p~ vertex displaced from the primary collision vertex. The displacement vec- 
tor between the p + p~ vertex and the primary vertex r is measured in the plane transverse 
to the beam direction. The most probable transverse b-hadron decay length in the laboratory 
frame H50j|5ll is calculated as 

L (4) 

where u is the unit vector in the direction of the ]/tp meson p~j and S _1 is the inverse of the 
sum of the primary and secondary vertex covariance matrices. From L xlJ the pseudo-proper 
decay length lyy = L xlJ my^l pj is computed as an estimate of the b-hadron decay length. The 
pseudo-proper decay length is measured with a resolution of ~ 35 pm. 

To measure the fraction of non-prompt ]/tp, the invariant-mass spectrum of p + p~ pairs and 
their ly^ distribution are fitted simultaneously using a two-dimensional unbinned maximum- 
likelihood fit in bins of pj, rapidity, and centrality with the fraction of non-prompt J/tp as a 
free parameter. The fitting procedure is similar to the one used in the pp analysis at y/s = 
7TeV [26] . The differences are: (i) the parametrisation of the lyy resolution function and (ii) the 
MC template used for the true tyy distribution of generated non-prompt ]/xp for which both 
muons have been reconstructed. Regarding (i), the reconstructed iyy distribution of simulated 
prompt ]/ip is better parametrised with a resolution function that is the sum of four Gaussians 
(the pp analysis at 7TeV used the sum of three Gaussians). Four of the eight fit parameters are 
fixed to the MC fit result and only the common mean, two widths, and one relative fraction 
are left free in the fits to the data. Regarding (ii), the lyy distribution of non-prompt ]/ip differs 
from that of the pp analysis because of the different heavy-ion tracking algorithm. In order 
to cope with the much higher detector occupancy, the PbPb tracking algorithm is done in one 
iteration and requires a pixel triplet seed to point to the reconstructed primary vertex within 
1 mm. Furthermore, the algorithm includes a filter at the last step that requires the track to 
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point back to the primary vertex within six times the primary vertex resolution. This reduces 
the reconstruction efficiency for J/ip with large values of £j/p i.e. it causes a difference in the 
prompt and non-prompt J/t/> reconstruction efficiencies that increases with the J/tp meson pj. 

The prompt J/i/> result is presented (in Section^ in the centrality bins 0-10%, 10-20%, 20-30%, 
30-40%, 40-50%, and 50-100%, while the non-prompt ]/xp result, given the smaller sample, 
is presented (in Section |7.2| in only two centrality bins, 0-20% and 20-100%. Examples of 
m^ l +y l - and lyy distributions are shown in Fig. |4] The two-dimensional fit results are shown 
as projections onto the mass and iyy axes. Integrated over centrality, the numbers of prompt 
and non-prompt ]/ip mesons with |y| < 2.4 and 6.5 < pj < 30 GeV/c are 307 ± 22 and 90 ± 13, 
respectively. 
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Figure 4: Invariant-mass spectra (left) and pseudo-proper decay length distributions (right) of 
pairs integrated over centrality (top) and for the 0-10% centrality bin (bottom). The spec- 
tra are integrated over the rapidity range |y| < 2.4 and the pj range 6.5 < pi < 30GeV/c. The 
projections of the two-dimensional fit onto the respective axes are overlaid as solid black lines. 
The dashed red lines show the fitted contribution of non-prompt ]/xp. The fitted background 
contributions are shown as dotted blue lines. 



In order to determine the systematic uncertainty on the yield extraction, the signal and back- 
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ground shapes are varied: for the signal mass shape, in addition to the default sum of the 
Crystal Ball and Gaussian functions, a single Gaussian and a single Crystal Ball function are 
tried. Alternatively, the oc and n parameters of the Crystal Ball function are fixed individu- 
ally for each pj and rapidity bin to the values found in the centrality integrated bin. This is 
in contrast to the default procedure in which the values for each rapidity bin are fixed to the 
values found in the bin integrated over centrality and all pj. For the background mass shape, a 
straight line is tried as an alternative. A crosscheck using a simple counting of the yield in the 
signal region after the subtraction of the same-sign spectrum leads to consistent results. The 
uncertainty on the fraction of non-prompt ]/ip due to the parametrisation of the iyy distribu- 
tion is estimated by varying the number of free parameters in the resolution function while the 
other parameters are fixed to their MC values. The systematic uncertainty is taken as the RMS 
of the yields obtained from the different variations of the fit function. The systematic uncer- 
tainties vary between 0.5% and 5.7% for the prompt j/ip yield, while the non-prompt ]/ip yield 
has uncertainties up to the extreme case of 14% in the most forward rapidity (1.6 < \y\ < 2.4) 
and lowest pj (3 < pj < 30 GeV/c) bin. 

4.2 Y(1S) Analysis 

To extract the Y(1S) yield, an extended unbinned maximum-likelihood fit to the invari- 
ant mass spectrum between 7 and 14 GeV/ c 2 is performed, integrated over pj, rapidity, and 
centrality, as shown in the left panel of Fig. [5] The measured mass line shape of each Y state is 
parametrised by a Crystal Ball function. Since the three Y resonances partially overlap in the 
measured dimuon mass spectrum, they are fitted simultaneously. Therefore, the probability 
distribution function describing the signal consists of three Crystal Ball functions. In addition 
to the three Y(nS) yields, the Y(1S) mass is the only parameter left free, to accommodate a pos- 
sible bias in the momentum scale calibration. The mass ratios between the states are fixed to 
their world average values fl48| , and the mass resolution is forced to scale linearly with the reso- 
nance mass. The Y(1S) resolution is fixed to the value found in the simulation, 92 MeV/c 2 . This 
value is consistent with what is measured when leaving this parameter free in a fit to the data, 
(122 ± 30) MeV/c 2 . The low-side tail parameters in the Crystal Ball function are also fixed to the 
values obtained from simulation. Finally, a second-order polynomial is chosen to describe the 
background in the mass range 7-14 GeV/c 2 . From this fit, before accounting for acceptance and 
efficiencies, the measured Y(1S) raw yield is 86 ± 12. The observed suppression of the excited 
states was discussed in |24|. The fitted mean value is mo = 

(9.441 ± 0.016) GeV/c 2 , which, for 
the same reason as for the J/i/7, is slightly below the PDG value m Y (is) = 9.460 GeV/c 2 lHH . 

The data are binned in pj and rapidity of the pairs, as well as in bins of the event cen- 

trality (0-10%, 10-20%, and 20-100%). The bins in rapidity are |y| < 1.2 and 1.2 < \y\ < 2.4. 
In contrast to the ]/xp case, CMS has acceptance for Y down to pj = OGeV/c over the full ra- 
pidity range. The pj bins in this analysis are < pj < 6.5GeV/c, 6.5 < pj < lOGeV/c, and 
10 < pj < 20GeV/c. There are only two events with a pair in the Y mass region and 

pj > 20 GeV/c. The invariant-mass distribution for the centrality bin 0-10% is illustrated in the 
right panel of Fig. [5] 

The systematic uncertainties are computed by varying the line shape in the following ways: 
(i) the Crystal Ball function tail parameters are varied randomly according to their covariance 
matrix and within conservative values covering imperfect knowledge of the amount of de- 
tector material and final-state radiation in the underlying process; (ii) the width is varied by 
±5 MeV/ c 2 , a value motivated by the current understanding of the detector performance; (iii) 
the background shape is changed from quadratic to linear, and the mass range of the fit is var- 
ied from 6-15 to 8-12GeV/c 2 ; the observed RMS of the results in each category is taken as the 
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CMS PbPbx/s^ 2.76 TeV 
Cent. 0-10%, |y| < 2.4 
< p T < 20 GeV/c 



92 MeV/c z (fixed to MC) - 




nv (GeV/c 2 



m„„ (GeV/cT 



Figure 5: Invariant-mass spectrum of }i + }i pairs (black circles) with pj < 20GeV/c and \y\ < 
2.4, for muons above 4GeV/c / integrated over centrality (left) and for the 0-10% centrality bin 
(right). 



systematic uncertainty. The quadratic sum of these three systematic uncertainties is dominated 
by the variation of the resolution of the mass fit, and is of the order of 10%, reaching 13% for 
the 0-10% centrality bin. As was the case for the ]/xp selection, a simple counting of the yield in 
the signal region after the subtraction of the same-sign spectrum leads to consistent results. 



5 Acceptance and Efficiency 
5.1 Acceptance 

The dimuon acceptance, A, is defined as the fraction of pairs for which both muons are 

declared detectable in the CMS detector with respect to all muon pairs produced in \y\ < 2.4, 

A( PT ,y;A e ) = ^^^ , (5) 

where: 

* ^detectable * s * ne numDer °f generated events in a given quarkonium (pj, y) bin in 
the MC simulation, for which both muons are detectable according to the selections 
defined in Eqs. ([T]) and (B; 

* ^generated is the number of all }i + pi~ pairs generated within the considered (pj, y) 
bin. 

The acceptance depends on the pj and y of the pair, and the polarization parameter Ag. 

Different polarizations of the ]/tp and Y(1S) will cause different single-muon angular distri- 
butions in the laboratory frame and, hence, different probabilities for the muons to fall inside 
the CMS detector acceptance. Since the quarkonium polarization has not been measured in 
heavy-ion or pp collisions at ^/s^ = 2.76 TeV, the prompt ]/tp and Y(1S) results are quoted for 
the unpolarized scenario only. For non-prompt J/xp the results are reported for the polariza- 
tion predicted by EVTGEN. The impact of the polarization on the acceptance is studied for the 
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most extreme polarization scenarios in the Collins-Soper and helicity frames. For fully longi- 
tudinal (transverse) polarized ]/xp in the Collins-Soper frame, the effect is found to be at most 
—20% (6%). In the helicity frame, the effects are at most 40% and —20% for the two scenar- 
ios. For Y (IS) the polarization effects range between —20% for longitudinal polarization in the 
Collins-Soper frame to 40% for transverse polarization in the helicity frame. 



The acceptance is calculated using the MC sample described in Section 3.1 The pj and rapidity 
dependencies of the ]/xp and Y(1S) acceptances are shown in Fig. [6] 
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Figure 6: Dimuon acceptance as a function of pj (left) and \y\ (right) for j/ip (red squares) 
and Y(1S) (green diamonds). Also shown in the right panel is the acceptance for }/ip with 
pj > 6.5 GeV/c (open black squares). The error bars represent the statistical uncertainties only. 

Since the acceptance is a function of both pj and y, uncertainties in the predicted distributions 
for these variables can lead to a systematic uncertainty in the average acceptance over a pj or 
y bin. To estimate these uncertainties, the shapes of the generated MC p T and y distributions 
are varied linearly by ±30% over the range \y\ < 2.4 and < pi < 30GeV/c (20GeV/c) for 
]/ip (Y(1S)). The RMS of the resulting changes in the acceptance for each pj and y bin are 
summed in quadrature to compute the overall systematic uncertainty from this source. The 
largest relative systematic uncertainties obtained are 4.2%, 3.2%, and 2.8% for the prompt ]/ip, 
non-prompt ]/tp, and Y(1S) acceptances, respectively. 

5.2 Efficiency 

The trigger, reconstruction, and selection efficiencies of pairs are evaluated using simu- 

lated MC signal events embedded in simulated PbPb events, as described in Section |3.1| The 
overall efficiency is calculated, in each analysis bin, as the fraction of generated events (pass- 
ing the single muon phase space cuts) where both muons are reconstructed, fulfil the quality 
selection criteria and pass the trigger requirements. In the embedded sample, the signal over 
background ratio is by construction higher than in data, so the background contribution un- 
derneath the resonance peak is negligible and the signal is extracted by simply counting the 
pairs in the quarkonium mass region. The counting method is crosschecked by using 
exactly the same fitting procedure as if the MC events were collision data. Only muons in the 
kinematic region defined by Eqs. ([I) and §2§ are considered. 

In Fig. [7j the efficiencies are shown as a function of the pair pj, y, and the event cen- 

trality, for each signal: red squares for prompt J/i/>, orange stars for non-prompt ]/tp, and green 
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5 Acceptance and Efficiency 



diamonds for Y(1S). As discussed in Section 4.1.2 the efficiency of non-prompt ]/rp is lower 
than that of prompt ]/ip, reaching about 35% for pj > 12GeV/c. The prompt j/ip efficiency 
increases with pj until reaching a plateau slightly above 50% at pj of about 12GeV/c, while 
the Y(1S) efficiency is ~55%, independent of pj. The efficiencies decrease slowly as a function 
of centrality because of the increasing occupancy in the silicon tracker; the relative difference 
between peripheral and central collisions is 17% for J/tp and 10% for Y(1S). The integrated 
efficiency values are 38.3%, 29.2%, and 54.5% for the prompt J/tp, non-prompt ]/xp (both with 
6.5 < pj < 30GeV/c, |y| < 2.4, and 0-100% centrality), and Y(1S) (with < p T < 20GeV/c, 
|y| < 2.4, and 0-100% centrality), respectively. 
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Figure 7: Combined trigger, reconstruction, and selection efficiencies as a function of quarko- 
nium pj and \y\, and event centrality, for each signal: red squares and orange stars for prompt 
and non-prompt ]/ip, respectively, and green diamonds for Y(1S). For better visibility, the 
prompt J/tp points are shifted by Apj = 0.5GeV/c, Ay = 0.05, and AN part = 2. Statistical (sys- 
tematic) uncertainties are shown as bars (boxes). The systematic uncertainties are the quadratic 
sum of the uncertainty on the kinematic distributions and the MC validation uncertainty. 

The systematic uncertainty on the final corrections due to the kinematic distributions is esti- 
mated by a ±30% variation of the slopes of the generated pj and rapidity shapes, similar to 
the acceptance variation described in the previous section. The systematic uncertainties are 
in the ranges 1.8-3.4%, 2.2-4.2%, and 1.4-2.7% for prompt J/xp, non-prompt ]/xp, and Y(1S), 
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respectively, including the statistical precision of the MC samples. 

The individual components of the MC efficiency are crosschecked using muons from ]/ip decays 
in simulated and collision data with a technique called tag-and-probe, similar to the one used 
for the corresponding pp measurement [26 j. In this method, high quality muons (the tags) 
are combined with muons that are selected without applying the selections whose efficiency 
is to be measured (the probes). Probe muons that fulfil these selections are then categorized 
as passing probes, the others as failing probes. A simultaneous fit of the two resulting invariant 
mass spectra (passing and failing) provides the efficiency of the probed selection. Because 
of correlations in the efficiency of matching silicon-tracker tracks to standalone muons, the 
total efficiency does not fully factorize into the individual components probed by this method. 
Therefore, the reconstruction and trigger efficiencies for }i + u~ pairs are directly obtained from 
the MC simulation, rather than as a product of the partial components. 

The fits are performed for tag-probe pairs with a pj above 6.5GeV/c as this is the region mea- 
sured over the full rapidity range, with and without applying the probed selection on one of 
the muons: 

1. The trigger efficiency is estimated by measuring the fraction of global muons (used as 
probes) associated to the double-muon trigger in an event sample selected by tag-muons 
associated to a single-muon trigger. A Crystal Ball function is used to describe the ]/xp 
peak. The p\ and rf l dependencies of the trigger efficiency are compatible between data 
and MC. For ]/xp with pj > 6.5 GeV/c, the p T and integrated trigger efficiency is 95.9% 
in MC and (95.1 ± 0.9)% in data. 

2. Standalone muons passing the quality selections required in this analysis are used to 
evaluate the efficiency of the silicon tracker reconstruction, which includes losses induced 
by the matching between the silicon-tracker track and the muon detector track, and by 
the imposed quality selection criteria (both on the global track and on its silicon-tracker 
segment). For this efficiency measurement, the signal is fitted with a Gaussian function 
and the background with a second-order polynomial. A Gaussian, rather than a Crystal 
Ball function, is used because of the poor momentum resolution of the standalone muons. 
No pj > 6.5GeV/c requirement was used, since the poorer momentum resolution of 
standalone muons would have biased the measurement. The single-muon efficiencies 
measured in MC and data of 84.9% and (83. 7^3)%, respectively, are in good agreement. 

The systematic uncertainty of the muon pair efficiency, 13.7%, is determined by comparing the 
tag-and-probe efficiencies evaluated in PbPb data and MC samples, and is dominated by the 
statistical uncertainties of the measurements. The standalone muon reconstruction efficiency 
(99% in the plateau) cannot be probed with silicon-tracker tracks because of the large charged 
particle multiplicity in PbPb collisions. Since this part of the reconstruction is identical to that 
used for pp data, a systematic uncertainty of 1%, reported in Ref. Il52l , is assumed. 

6 The pp Baseline Measurement 

A pp run at y/s = 2.76 TeV was taken in March 2011. The integrated luminosity was 231 nb _1 , 
with an associated uncertainty of 6%. For hard-scattering processes, the integrated luminosity 
of the pp sample is comparable to that of the PbPb sample (7.28 f/b -1 • 208 2 « 315 nb _1 ). 

Given the higher instantaneous luminosity, the Level-1 trigger required slightly higher quality 
muons in the pp run than in the PbPb run. The offline event selection is the same as in the PbPb 
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6 The pp Baseline Measurement 



analysis, only slightly relaxed for the HF coincidence requirement: instead of three towers, only 
one tower with at least 3 GeV deposited is required in the pp case. The same reconstruction 
algorithm, i.e. the one optimized for the heavy-ion environment, is used for both pp and PbPb 
data. The products of the trigger, reconstruction, and selection efficiencies determined in pp 
MC simulations are 42.5%, 34.5%, and 55.1% for the prompt ]/ip, non-prompt ]/ip (both with 
6.5 < p T < 30GeV/c, |y| < 2.4), and Y(1S) (withO < p T < 20GeV/c, |y| < 2.4), respectively. 

The accuracy of the MC simulation in describing the trigger efficiency is crosschecked with the 
tag-and-probe method in the same way as for the PbPb analysis discussed in Section 5.2 For 
muons from decays of J/tp with pj > 6.5 GeV/c, the p l j and integrated trigger efficiencies are 
(92.5 ± 0.6)% in data and (94.3 ± 0.2)% in MC. In the same phase-space, the tracking and muon 
selection efficiency is (82.5 ± 2.4)% in data and (84.6 ± 1.0)% in MC. For the standalone muon 
reconstruction efficiency a systematic uncertainty of 1% is assigned, as reported in Ref. [52]. As 
in the PbPb case, the systematic uncertainty of the muon pair efficiency in pp collisions, 13.7%, 
is determined by comparing the tag-and-probe efficiencies evaluated in data and MC samples, 
and is dominated by the statistical uncertainties of the measurements. 

The quarkonium signals in pp collisions are extracted following the same methods as in PbPb 



collisions, described in Sections 4.1 and 4.2 apart from the non-prompt ]/xp signal extraction: 
the four Gaussians of the lifetime resolution are fixed to the MC values because of the lack of 
events in the dimuon mass sidebands. The systematic uncertainty on the signal extraction in 
pp is 10% for Y(1S) and varies, depending on pj and rapidity, between 0.4 and 6.2% for prompt 
]/ip and between 5 and 20% for non-prompt ]/tp. The fit results for the prompt and non-prompt 
]/ip yield extraction are shown in Fig. [8]for \y\ < 2.4 and 6.5 < pj < 30GeV/c. The numbers of 
prompt and non-prompt ]/ip mesons in this rapidity and pj range are 820 ± 34 and 206 ± 20, 
respectively. 
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Figure 8: Non-prompt J/t/> signal extraction for pp collisions at \fs = 2.76 TeV: dimuon invari- 
ant mass fit (left) and pseudo-proper decay length fit (right). 



The invariant-mass spectrum of pairs in the Y region from pp collisions is shown in 

Fig. |9j The same procedure as the one described for the PbPb analysis is used. The number of 
Y(1S) mesons with \y\ < 2.4 and < p T < 20GeV/c is 101 ± 12. The fit result of the excited 
states is discussed in [24] . 

The differential cross section results include the systematic uncertainties of the reconstruction 
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Figure 9: The pp dimuon invariant-mass distribution in the range pj < 20GeV/c for 
and the result of the fit to the Y resonances. 
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efficiency and acceptance, estimated in the same way as for the PbPb analysis. The systematic 
uncertainties on the efficiencies are 1.6-3%, 1.4-2%, and 0.4-0.9% for prompt ]/ip, non-prompt 
]/ip, and Y(1S), respectively. The uncertainty on the acceptance is identical in the pp and PbPb 
analyses. 

For the measurement of the nuclear modification factors, in which the ratio of PbPb to pp 
results is computed, most of the reconstruction systematic uncertainties cancel out because the 
same algorithm is used. However, the following factors must be accounted for: 

1. The luminosity uncertainty. This is a global systematic uncertainty of 6% that allows all 
measured nuclear modification factors to change by a common scale-factor. Since the 
PbPb yield is normalized by the number of minimum-bias events, which has a negligible 
uncertainty, no systematic uncertainty on the PbPb luminosity has to be considered. 

2. The uncertainty on Taa- F° r results integrated over centrality this is a global systematic 
uncertainty of 5.7%, based on the Glauber model employed. For results as a function of 
centrality, the uncertainty varies between a minimum of 4.3% in the most central bin and 
a maximum of 15% in the most peripheral bin [36J. 

3. The systematic uncertainty associated with the trigger efficiency. The ratios between the 
tag-and-probe efficiencies obtained in pp and PbPb are the same in data and MC events, 
within the statistical accuracy of the data (1% for the single-muon efficiency). Twice this 
value (2%) is assigned as the uncertainty on the difference of the trigger efficiencies of 

pairs in PbPb and pp collisions. 

4. The tracking efficiency uncertainty due to different charged particle multiplicities in pp 
and PbPb collisions. The ratios between the tag-and-probe efficiencies obtained in pp and 
central PbPb events are the same in data and MC events, within the statistical accuracy 
of the data (6.8% for the single-muon efficiency). This value is propagated as the tracking 
systematic uncertainty in all the ratios of PbPb to pp data. 
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7 Results 



7 Results 

The double-differential quarkonium cross sections in PbPb collisions are reported in the form 

1 d 2 N 1 1 N m 



T AA dydp T T aa N M b Ai/Ap T Ae 

while in pp collisions they are calculated as 

d 2 a 1 1 N Q q 
di/dpx £ pp AyAp T Ae ' 

where: 



(6) 



(7) 



• ^qq * s me number of measured prompt ]/ip, non-prompt ]/tp, or Y(1S) in the }i + }i 
decay channel; 

• Nmb is the number of minimum-bias events sampled by the event selection; when 
binned in centrality only the fraction of minimum-bias events in that centrality bin 
is considered; 

• A is the geometric acceptance, which depends on the pj and y of the quarkonium 
state; 

• £ is the combined trigger and reconstruction efficiency, which depends on the pj and 
y of the quarkonium state and on the centrality of the collision; 

• Ay and Apj are the bin widths in rapidity and pj, respectively; 

• T AA is the nuclear overlap function, which depends on the collision centrality; 

• Cpp = (231 ± 14) nb _1 is the integrated luminosity of the pp data set. 

Following Eq. the uncorrected yields of inclusive, prompt and non-prompt ]/ip, and Y(1S), 
measured in PbPb collisions are corrected for acceptance and efficiency, and converted into 
yields divided by the nuclear overlap function T AA . These quantities can be directly compared 
to cross sections in pp collisions measured from the raw yields according to Eq. 0. The rapidity 
and centrality-dependent results are presented integrated over pj. All results are presented for 
the unpolarized scenario only and are tabulated in Tables |4]-[7] of Appendix |A| 

The systematic uncertainties detailed in the previous sections are summarized in Tables [2]and [3] 
The relative uncertainties for all terms appearing in Eqs. (|6j and |7|) are added in quadrature, 
leading to a total of 15-21% on the corrected yields. For results plotted as a function of pj 
or rapidity, the systematic uncertainty on T AA enters as a global uncertainty on the scale and 
is not included in the systematic uncertainties of the yields. As a function of centrality, the 
uncertainty on T AA varies point-to-point and is included in the systematic uncertainties of the 
yields. 

The nuclear modification factor, 

7? _ PP 

Np b p b (QQ) 

PP CQ\ 

T aa N M b Npp(QQ) epbPb 

is calculated from the raw yields Npbpb(QQ) and N pp (QQ), correcting only for the multiplicity- 
dependent fraction of the efficiency (^j^ ~ 1-16 for the most central bin); the pj and rapid- 
ity dependencies of the efficiency cancel in the ratio. These results are also tabulated in Ap- 
pendix [A] 



7.1 Inclusive and Prompt J/ip 
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Table 2: Point-to-point systematic uncertainties on the prompt ]/ip, non-prompt J/t/;, and Y(1S) 
yields measured in PbPb collisions. 





prompt J/t/; (%) 


non-prompt J/t/; (%) 


Y(1S) (%) 


Yield extraction 


0.5-5.7 


1.5-14.0 


8.7-13.4 


Efficiency 


1.8-3.4 


2.2-4.2 


1.4-2.7 


Acceptance 


0.9-4.2 


2.0-3.2 


1.5-2.8 


MC Validation 


13.7 


13.7 


13.7 


Stand-alone ]i reco. 


1.0 


1.0 


1.0 


Taa 


4.3-15.0 


4.6-8.6 


4.3-8.6 



Total 15-21 15-21 18-20 



Table 3: Point-to-point systematic uncertainties on the prompt J/t/;, non-prompt J/t/;, and Y(1S) 
yields measured in pp collisions. 





prompt J/t/; (%) 


non-prompt J/t/; (%) 


Y(1S) (%) 


Yield extraction 


0.8-5.3 


5.3-16.8 


10.0 


Efficiency 


1.6-3.0 


1.4-2.0 


0.4-0.9 


Acceptance 


0.9-4.2 


2.0-3.2 


1.5-2.8 


MC Validation 


13.7 


13.7 


13.7 


Stand-alone \i reco. 


1.0 


1.0 


1.0 


Total 


14-16 


15-22 


17-18 



In all figures showing results, statistical uncertainties are represented by error bars and system- 
atic uncertainties by boxes. Results as a function of rapidity are averaged over the positive and 
negative rapidity regions. 

7.1 Inclusive and Prompt J/ip 

The inclusive and prompt J/t/; differential yields in PbPb collisions, divided by Taa, are shown 



in the left panel of Fig. 10 as a function of pj, for |y| < 2.4 and integrated over centrality. The 
corresponding pp cross sections are also shown. The suppression of the prompt J/t/; yield by 
a factor of ~ 3 with respect to pp is easier to appreciate through the Raa observable, shown 



in the right panel of Fig. 10 The Raa measurements do not exhibit a pj dependence over the 



measured pj range, while there is an indication of less suppression in the most forward rapidity 



bin (1.6 < \y\ < 2.4) in comparison to the mid-rapidity bin, as shown in Fig. 11 At forward 
rapidity, in addition to 6.5 < pj < 30GeV/c the nuclear modification factor is measured for 
lower pj (down to 3 GeV/c) without observing a significant change, as can be seen in Table |5j 

The inclusive J/t/; yield in PbPb collisions divided by Taa, integrated over the pj range 6.5- 



30 GeV/c and \y\ < 2.4, is shown in the left panel of Fig. 12 as a function of N part . Also included 
is the prompt J/t/; yield, which exhibits the same centrality dependence as the inclusive J/t/;: 
from the 50-100% centrality bin ((N part ) = 22.1) to the 10% most central collisions ((N part ) = 
355.4) the yield divided by Taa falls by a factor of ~ 2.6. The results are compared to the cross 
sections measured in pp, showing that prompt J/t/; are already suppressed in peripheral PbPb 
collisions. The Raa of prompt J/t/; as a function of N part is shown in the right panel of Fig. 12 a 
suppression of ~ 5 is observed in the 10% most central PbPb collisions with respect to pp. This 
suppression is reduced in more peripheral collisions, reaching a factor of ~ 1.6 in the 50-100% 
centrality bin. 
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7.2 Non-prompt }lxp 

The uncorrected fraction of non-prompt }/\p is obtained from the two-dimensional fit to the 
invariant mass and lyq spectra discussed in Section 4.1.2 To obtain the corrected b fraction, 
which is the ratio of non-prompt to inclusive J/i/?, the raw fraction is corrected for the different 
reconstruction efficiencies and acceptances for prompt and non-prompt }/ip. The b fraction 
in pp and in PbPb (integrated over centrality) at ^s NN = 2.76 TeV is presented in Fig. Il3| as a 
function of pj, for several rapidity bins, together with results from CDF [41] and CMS 11261 at 
other collision energies. There is good agreement, within uncertainties, between the earlier 
results and the present measurements. 

The non-prompt ]/xp yield in PbPb collisions divided by Taa, integrated over the pj range 6.5- 
30GeV/c and \y\ < 2.4, is shown in the left panel of Fig. 14 as a function of Np ar t, together 
with the pp cross section. Non-prompt }/ip are suppressed by a factor of ~ 2.6 with respect to 



pp collisions, as can be seen in the right panel of Fig. 14 The suppression does not exhibit a 
centrality dependence, but the most peripheral centrality bin (20-100%, (Np art ) = 64.2) is very 
broad. Hard processes, such as quarkonium and b-hadron production, are produced following 
a scaling with the number of nucleon-nucleon collisions, thus most events in such a large bin 
occur towards its most central edge. 

7.3 Y(1S) 



In Fig. 15 the Y(1S) yield divided by Taa i n PbPb collisions and its cross section in pp collisions 
are shown as a function of pj; the Raa of Y ( IS) is displayed in the right panel of Fig. 15 The pj 
dependence shows a significant suppression, by a factor of ~ 2.3 at low pj, that disappears for 
pj > 6.5GeV/c. The rapidity dependence indicates a slightly smaller suppression at forward 



rapidity, as shown in Fig. 16 However, the statistical uncertainties are too large to draw strong 



20 



7 Results 



c 1 
o 

o 0.9 

CO 

S 0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



Q 



+ 

A 







1 I 1 1 1 1 I 
CMS PbPb 

CMS PbPb 
CMS pp Vs 
CMS pp 
CMS pp 
CMS pp 
CMS pp 
CDF pp 



T 



T 



T 



s M =2.76TeV(|y|<2.4) 
s NN = 2.76TeV(1.6<|y|<2.4) 
= 2.76 TeV (|y|<2.4) 
= 2.76 TeV (1.6<|y|<2.4) 
= 7TeV(1.6<|y|<2.4) 
= 7TeV (1.2 < |y| < 1.6) 
= 7TeV(|y| < 1.2) 
1.96 TeV (|y[<0.6) 



10 



15 20 

(GeV/c) 



25 



30 



Figure 13: b fraction of ]/ip production in pp and PbPb collisions at ^s NN = 2.76 TeV as a 
function of pj for the rapidity bins \y\ < 2.4 and 1.6 < |y| < 2.4, compared to b fractions 
measured by CDF in pp collisions at \/s = 1.96 TeV [41J and by CMS in pp collisions at a/s = 
7 TeV [26 1. Points are plotted at their measured average pj. Statistical (systematic) uncertainties 
are shown as bars (boxes). 



o 
z 

"D 



2.5 



1.5 



0.5 





1 1 - 


- CMS pp & PbPb \ls m 


= 2.76 TeV - 


★ Non-prompt JA|/ (PbPb) 


+ (PP) 


i 

1 20-1 00% 


0-20% Z 


- i 

I lyl<2.4 


+ J 


- 6.5 < p T < 30 GeV/c 




I , , , , I , , , , 1 , , , , 1 , , , , I , , , , I , 


, , 1 , , , , 1 , , , f 



50 100 150 200 250 300 350 400 



1.2 



0.6 



0.4 



7 CMS PbPb \ls m 

★ Non-prompt JA|f 


= 2.76 TeV " 


: 20-100% 


0-20% : 


\ \ 

r lyl < 2.4 


+ i 


- 6.5 <p <30GeV/c 

"t 1 T t I T t T T I T T t T I t T T t I 1 





50 100 150 200 250 300 350 400 



N, 



part 



Figure 14: Left: non-prompt ]/ip yield divided by (orange stars) as a function of N par t 
compared to the non-prompt J/ip cross section measured in pp (black cross). Right: nuclear 
modification factor Raa of non-prompt ]/xp as a function of N par t- A global uncertainty of 6%, 
from the integrated luminosity of the pp data sample, is shown as a grey box at Raa = 1- 
Statistical (systematic) uncertainties are shown as bars (boxes). 



21 



conclusions on any pj or rapidity dependence. The Y(1S) yield in PbPb collisions divided 
by Taa and the Y(1S) Raa are presented as a function of N par t in the left and right panels of 
Fig. 17 respectively. Within uncertainties, no centrality dependence of the Y(1S) suppression 
is observed. 
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8 Discussion 

This paper has presented the first measurements of the prompt and non-prompt ]/ip, as well 
as the Y(1S) mesons, via their decays into pairs in PbPb and pp collisions at ^/s NN = 

2.76 TeV. The results are based on data recorded with the CMS detector from the first LHC 
PbPb run in 2010, and from a pp run during March 2011 at a/s = 2.76 TeV. 

The prompt ]/\p cross section shows a factor of two suppression in central PbPb collisions with 
respect to peripheral collisions for ]/ip with 6.5 < pj < 30GeV/c. With respect to pp, a nuclear 
modification factor of Raa = 0.20 ± 0.03 (stat.) ± 0.01 (syst.) has been measured in the 10% 
most central collisions. Prompt ]/tp produced in peripheral collisions are already suppressed 
with respect to pp: Raa = 0.61 ± 0.12 (stat.) ± 0.10 (syst.) in the 50-100% centrality bin. While 
no pj dependence is observed in the measured pj range, within uncertainties, less suppression 
is observed at forward rapidity (Raa = 0.43 ± 0.06 (stat.) ± 0.01 (syst.)) than at mid-rapidity 
(R AA = 0.29 ± 0.04 (stat.) ± 0.02 (syst.)). 

A comparison of the Raa centrality dependence to results measured for pj < 5GeV/c by 
PHENIX [21J in AuAu collisions at ^/s NN = 200 GeV shows a similar suppression, despite the 
different collision energies and kinematic ranges. Integrated over centrality, CMS has mea- 
sured an inclusive }/ip nuclear modification factor of Raa = 0.41 ± 0.05 (stat.) ± 0.02 (syst.) in 
the most forward rapidity bin (1.6 < \y\ < 2.4) in the pj range 3 < pj < 30GeV/c. 



22 



8 Discussion 



>, 1.6 
"a 

-a 1.4 



1.2 

1 

0.8 
0.6 
0.4 
0.2 



L CMS pp & PbPb 


\/s NN = 2.76TeVi 


♦ T(1S)(PbPb) 

i — i — - 


+ (PP) 

-*-4 


: Cent. 0-100% 

-, M I , , I I I , , I I , , I I I , I I I , I M 


< p < 20 GeV/c 

1 1 T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



CC 1.4 



1.2 



0.8 
0.6 
0.4 
0.2 



7 CMS PbPb * 


^s NN = 2.76TeV ~ 


♦ T(1S) 














i 1 








> 


- Cent. 0-100% 

- 0<p <20GeV/c 

' t I , , , I , , , I , i , I , , 1 1 


— i 


, i , , , i , , , i , , , i , , r 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 



|y| |y| 
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The global scale uncertainties on the PbPb data due to Taa (5.7%) and the pp integrated lumi- 
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Figure 17: Left: Y(1S) yield divided by Taa (green diamonds) as a function of N par t compared 
to the Y(1S) cross section measured in pp (black cross). Right: nuclear modification factor Raa 
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A strong suppression of non-prompt J/tp mesons is observed in PbPb collisions when compared 
to pp collisions. This is the first unambiguous measurement of b-hadron suppression in heavy- 
ion collisions, which is likely connected to in-medium energy loss of b quarks. The average pj 
of the non-prompt J/xp in the measured kinematic range is ~ 10 GeV/c. Based on simulations of 
b-hadron decays, this translates into an average b-hadron pj of ~ 13GeV/c. The suppression 
of non-prompt ]/ip is of a comparable magnitude to the charged hadron Raa measured by AL- 
ICE H53H . which reflects the in-medium energy loss of light quarks. The non-prompt ]/xp yield, 
though strongly suppressed (Raa = 0-37 ± 0.08 (stat.) ± 0.02 (syst.)) in the 20% most central 
collisions, shows no strong centrality dependence, within uncertainties, when compared to a 
broad peripheral region (20-100%). Furthermore, this suppression of non-prompt ]/tp is com- 
parable in size to that observed for high-px single electrons from semileptonic heavy-flavour 
decays at RHIC II28TI301 in which charm and bottom decays were not separated. 

The Y(1S) yield divided by Taa as a function of px, rapidity, and centrality has been measured 
in PbPb collisions. No strong centrality dependence is observed within the uncertainties. The 
nuclear modification factor integrated over centrality is Raa = 0-63 ± 0.11 (stat.) ± 0.09 (syst.). 
This suppression is observed predominantly at low pj. Using pp collisions at a/ s = 1-8 TeV, 
CDF measured the fraction of directly produced Y(1S) as (50.9 ± 8.2 (stat.) ± 9.0 (syst.))% for 
Y(1S) with pj > 8GeV/c [54|. Therefore, the Y(1S) suppression presented in this paper could 
be indirectly caused by the suppression of excited Y states, as indicated by earlier results from 

cms m. 

9 Summary 

In summary, CMS has presented the first measurements of prompt ]/ip, non-prompt ]/rp, and 
Y(1S) suppression in PbPb collisions at ^/s NN = 2.76 TeV. Prompt ]/tp are found to be sup- 
pressed, with a strong centrality dependence. By measuring non-prompt ]/tp, CMS has directly 
observed the suppression of b hadrons for the first time. The measurement of Y(1S) suppres- 
sion, together with the suppression of the Y(2S+3S) states |24j, marks the first steps of detailed 
bottomonium studies in heavy-ion collisions. 
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A Tables of Results 

Table 4: Yield per unit of rapidity of inclusive J/tp divided by Taa an d nuclear modification 
factor Raa as a function of ]/ip rapidity, pj, and collision centrality. The average pj value for 
each bin is given. Listed uncertainties are statistical first, systematic second, and global scale 
third. The latter includes the uncertainties on Taa an d on the pp integrated luminosity. 
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Table 5: Yield per unit of rapidity of prompt J/rp divided by Taa and nuclear modification factor 
Raa a s a function of ]/ip rapidity pj, and collision centrality. The average pj value for each bin 
is given. Listed uncertainties are statistical first, systematic second, and global scale third. The 
latter includes the uncertainties on Taa an d on the pp integrated luminosity. 
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Table 6: Yield per unit of rapidity of non-prompt J/tp divided by Taa an d nuclear modification 
factor R A a a s a function of J/ip rapidity, pj, and collision centrality. The average pj value for 
each bin is given. Listed uncertainties are statistical first, systematic second, and global scale 
third. The latter includes the uncertainties on Taa and on the pp integrated luminosity. 
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Table 7: Yield per unit of rapidity of Y(1S) divided by Taa and nuclear modification factor Raa 
as a function of Y(1S) rapidity pj, and collision centrality. The average px value for each bin 
is given. Listed uncertainties are statistical first, systematic second, and global scale third. The 
latter includes the uncertainties on Taa and on the pp integrated luminosity. 
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